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ABSTRACT 


The Department of Aeronautics and Astronautics at the Naval Postgraduate School (NPS) is 
expanding its helicopter research capabilities in order to facilitate present and future research 
demands. The rapidly changing needs have already out paced available assets. Therefore it was 
necessary to design and develop a new remotely piloted helicopter (RPH) that would meet present 
needs, NOTAR and HHC, and be flexible enough to meet future needs. The research efforts 
encompassed by this thesis are defining the present needs, investigating what type/size of RPH would 
fulfill these needs, procuring this asset, and analyzing its capabilities. Based on a defined payload, 
helicopter trends are analyzed to determine an estimate of the overall RPH size (gross weight) and 
engine size required. A preliminary design process validates these figures. Choosing to procure an 
RPH instead of building one, a detailed performance analysis is conducted on the main rotor system. 
This analysis includes blade vibration analysis, retreating blade stall analysis, and power required 
analysis. Modification of the RPH’s main rotor hub, drive train, and landing gear are studied and 
recommendations presented. This research effort is a continuation of a long- term program to provide 


NPS with robust assets to support present and future rotorcraft research efforts. 


ill 


TABLE OF CONTENTS 


El INTRODUCTION 


18 Ee BACKGROUND 
A. HIGHER HARMONIC CONTROL (HHC) 
B. NOTAR 


C. AUTONOMOUS LANDING AND TAKEOFF SYSTEM (ALTOS) 


idole DEVELOPING SCALE MODEL TRENDS 


AISI PRELIMINARY DESIGN 
A. CONSTRAINT DIAGRAM 
Be PRELIMINARY ESTIMATES 
1. Gross Weight and Blade Radius 
2. One Hour Estimate 


) 


a. Maximum Forward Velocity (V,,, 
by TROPPO SazinG 
c. Staneing Permaer 


3. DESIGN FROM SCRATCH OR PROCURE 
Ve HUMMINGBIRD ACQUISITION 
A. CAPABILITIES 


B. SHORTFALLS 


ab 


14 


14 


bes. 


LS 


18 


Bl 


PAN, 


Pa AL 


Ze 


23 


2a 


25 


Wal 


pair. 


A. 


Wert 


1. Two Blade Rotor 
2. Lack of Autorotational Capability 


3. Low Reynolds Number 


BLADE DESIGN 


SPECIFICS 
1. AIRFOIL 
2. IWiSt 


32 PEP VEnOCELY 

TRADEOFF 

SL abil Ty 

BLADE DYNAMICS 

jeer COdUeE TON 

2. Main Rotor Blades 

3. Myklestad Determinant Method 
4. Blade Modes 


5. Southwell Plot 


PERFORMANCE ANALYSIS 


POWER CALCULATIONS 
1. Power Required 
2. Power Available 


RETREATING BLADE STALL 


MODIFICATIONS 


LANDING GEAR 


20 


Ze 


ae 


ag 


29 


29 


30 


ape 


a2 


35 


39 


32 


40 


45 


oy 


54 


60 


60 


60 


63 


64 


seal 


dak 


B. DRIVE TRAIN MODIFICATIONS 
Engine to Main Rotor Gear Ratio 
Freewheeling Unit 
AX. CONCLUSIONS AND RECOMMENDATIONS 
A. CONCLUSIONS 
B. RECOMMENDATIONS 
APPENDIX MYKLESTAD FORTRAN CODE 
APPENDIX MYKLESTAD NATURAL FREQUENCIES 
APPENDIX MODE SHAPE FORTRAN CODE 
APPENDIX PRELIMINARY POWER CALCULATIONS 
APPENDIX BLADE TIP PITCH ANGLE WORKSHEET 


LIST OF REFERENCES 


INL TIAL DISTRIBULICN Eros 


Wea 


wan 


gil 


TZ 


73 


ve 


TS 


77 


80 


88 


Se) 1. 


gz 


he) 


pe) 


I. INTRODUCTION 


Helicopter research in the Department of Aeronautics and 
Astronautics at the Naval Postgraduate School (NPS) is 
rapidly accelerating. This fast pace has created a need for 
expanding the capabilities of the department’s Remotely 
Piloted Vehicles (RPV). Valuable research has been 
accomplished through the employment of the existing RPVs. 
The GMP Legend, a commercially produced radio controlled 
(RC) helicopter, was used to produce baseline vibration 
analysis and to validate measurement techniques. The 
results from this work will be of great value in subsequent 
Higher Harmonic Control (HHC) research. The Bruiser, a 
limited production, 20 pound payload RPV developed by 
Pacific RPV, was used ina shake test to obtain an airframe 
modal analysis. 

The preceding research has been of great value to the 
department, but in order to advance beyond proof of concept 
and into scale model analysis, a new helicopter was 
required. Incorporated in the requirements were a 
helicopter with a tip speed and chord which produced a 
Reynolds number that was within an acceptable range for 
comparison with a full scale helicopter; a helicopter of a 


size large enough to carry a No-Tail-Rotor (NOTAR) tail boom 


that would also be testable in the quarter scale (14 by 22 
foot) wind tunnel at NASA Langley; and a helicopter with 
enough flexibility to be adaptable to future research needs. 
This impetus developed the challenge of designing a suitable 
RPV for the department. 

A critical amount of background knowledge and direction 
was obtained on a research trip to the Aerostructures 
Directorate at NASA Langley. They were heavily involved in 
the design, manufacture and test of their own Free-Flight 
Rotorcraft Research Vehicle (FFRRV), which was directly 
along the same lines of interest of NPS. Their RPV was 
twice the size and considerably more complex than that 
desired by NPS. The head of the directorate, Arthur E. 
Phelps III, was invaluable in passing on their corporate 


knowledge, saving untold hours in achieving the NPS goal. 


ig a BACKGROUND 


A. HIGHER HARMONIC CONTROL (HHC) 

The control of vibrations has always been of great 
concern to both the helicopter designer and to the 
helicopter pilot. It has been a continuing source of 
agitation throughout the years and a focus of enormous 
amounts of research assets. The current means by which the 
vibrations are reduced are through passive devices which 
either isolate the source of vibration (isolators) or 
diffuse the vibration level (absorbers). These vibration 
absorption mechanisms are usually restricted to a narrow 
scope of flight conditions and vibratory frequencies. The 
use of HHC, a relatively new technology, is an active 
vibration reduction device, vice the passive ones just 
mentioned. HHC functions by altering the aerodynamic loads 
on the rotor, and therefore the vibratory forces and moments 
which cause the airframe to vibrate are reduced. [Ref. 1] 

In earlier full scale HHC testing on the Hughes 
Helicopter OH-6A, it was determined that not only were the 
vibrations successfully reduced, but the performance of the 
helicopter was also improved. Other helicopter companies 
tried to duplicate this resulting improved performance but 


were unsuccessful. Corroborating this performance 


enhancement locally would have great merit. The RPVs on 
hand did not allow this type of research to be conducted 
because they could not produce Reynolds numbers over the 
Main rotor blade that would allow data comparison with the 
data obtained on the OH-6A. A helicopter with a tip speed 
and chord length that would produce Reynolds numbers on the 
order of two million was required. This shortcoming 


provided some of the drive to obtain a new helicopter RPV. 


B. NOTAR 

Presently, Lt. M. Borno, in conjunction with McDonnell 
Douglas Helicopter Corporation, is conducting thesis 
research which will produce a fully operational scale model 
NOTAR tailboom. The tailboom requirement was that it was to 
be big enough to test in the quarter scale wind tunnel at 
NASA Langley. The current RPVs lacked the size to 
mechanically support and drive a tailboom of this scale. 
This was a second driver behind acquiring a new helicopter 


Rev; 


C. AUTONOMOUS LANDING AND TAKEOFF SYSTEM (ALTOS) 

The United States Navy’s growing interest in RPVs and 
their fleet applications has now reached a point which 
requires significant background research. As the Navy leans 
toward fully autonomous RPVs, it lacks experience and 


expertise in the most effective and efficient means to 


launch and recover these vehicles. Recently, Orion Aviation 
entered into a contract with the United States Navy in 
conjunction with NPS to develop and test five potential 
ALTOS concepts and have them analyzed for merit by students 
at NPS. Once the best ALTOS concept is determined, it will 
be built by Orion Aviation and then will be demonstrated 
using an NPS RPV. The expanded capabilities that the new 
NPS RPV would add to the existing resources would provide 
greater flexibility in the demonstration and validation 


phase of the ALTOS program. 


Lit: DEVELOPING SCALE MODEL TRENDS 


There is a great deal of documentation for the design of 
full scale helicopters regardless of weight range, but there 
exits little or no documentation for design of a scale model 
helicopter. It was necessary to determine trends for gross 
weight, takeoff weight, disk loading (DL), rotor radius, 
solidity (0), and blade loading (BL) from full scale 
helicopters in order to determine what the values should be 
for quarter-scale size. The trends that were developed were 
typically linear with an adjustment of either the x- 
intercept or the y-intercept usually required. The 
adjustment was made assuming the trends held true for all 
weight ranges, but at the lower weights the line had to be 
shifted to accommodate scale sizes. The trends that were 
generated show this assumption to be good. 

The first trend determined was between the load a 
helicopter could carry versus the gross weight of the 
helicopter. Fig. 1 shows data for thirteen different full 
scale helicopters plus the Bruiser. 


A linear trend analysis produced Eqn. 1: 


Gross Weight (1) 


Payload = 270 + 
Ao 


The slope of the line was assumed to be true for all weight 


ranges, but the y-intercept was decreased for gross weight 
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Figure 1 Payload vs Gross Weight 


less than 500 pounds. The final equation that was 


determined is shown in Eqn. 2. 


Gross Weight = 4.52 + Payload (2) 


The needs as stated in the previous chapter required a 
payload of between 20 and 30 pounds. This equates to a 


helicopter between 90 and 135 pounds. 


The next trend that was studied was takeoff weight 
versus usable power. The desire was to use this trend to 
provide a general idea of what power would be required based 
on the above weight range. The trend shown in Fig. 2 was 


used to derive Eqn. 3. 


Takeoff Weight 


7.38 oa 


Usable Power = -70.66 + 
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Takeoff Wt vs. Usable Power 
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Usable Pwr = TO Wit/7 38 






Figure 2 Usable Power vs Takeoff Weight 


Again, assuming the slope was accurate over the entire 
weight range, the intercept point was dropped to zero. The 
final equation is shown below. This equation shows that 


power required would range between 12 to 18 horsepower. 


Takeoff Weight 


7.38 oy 


Usable Power = 


Next, the main rotor blade solidity and blade loading (Ct/c) 
versus takeoff weight were compared. Fig. 3 shows that main 
rotor blade solidity is independent of takeoff weight. The 
minimum solidity was found to be 0.03 for the Robinson R-22 
Beta. The maximum solidity was 0.098 for the Aerospatiale 
AS 332 Ll. The average solidity was approximately 0.07. 
Fig. 4 shows a similar trend to that of the blade solidity. 
The blade loading was also independent of takeoff weight. 
The minimum was 0.05 for the Bell 412. The maximum was 
Seo95 Cor 
the Aerospatiale SA 365N, and the average was 0.075. 

The typical disk loading (DL) for model helicopters is 
between 1.0 and 2.0. Fig. 5 shows that the model helicopter 
does not fall within the normal range of full scale 


helicopters. The minimum full scale DL of 2.8 is that of 
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Takeoff Wt vs. Blade Solidity 
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Solidity vs Takeoff Weight 
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Takeoff Wt vs. Blade Loading 
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Figure 4 Blade Loading vs Takeoff Weight 
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Figure 5 Disk Loading vs Takeoff Weight 


the Robinson R-22. Full scale trends were inconclusive in 
providing guidelines for the model helicopter. Fig. 6 shows 
how power loading (PL) varies with disk loading. Power 
loading tells how much weight can be lifted for a given 
horsepower. The lower the DL, the greater the PL that can 
be achieved. In summary, the trend analysis shows that a 
desired payload of 20 to 30 pounds will require a helicopter 


weight between 90 and 135 pounds. An engine of 12 to 18 
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Figure 6 Power Loading vs Disk Loading 


horsepower would be required for a helicopter in this weight 
range. The solidity should be in the vicinity of 0.07 and 
the BL should be in the vicinity of 0.075. These trends 
provide a good starting point for the design process anda 


very good point of comparison for the preliminary design 


results. 
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IV. PRELIMINARY DESIGN 


A. CONSTRAINT DIAGRAM 

The first task to accomplish was to construct a rotor 
blade tip speed constraint diagram. Based on industry 
criteria (Ref. 2:p. 90], a tip speed constraint diagram was 
made using the following limits. The upper boundary can 
either be set by a noise limit which is 750 fps or a hover 
tip mach number limit which is M<0.69 (771 fps). The hover 
tip mach number was used for this constraint diagram. The 
lower limit is set by requirement to store kinetic energy in 
the rotor system in case of power failure, in other words an 
autorotational limit. This limit is set at 400 fps. Two 
other limits that are important to consider are that of the 
advancing blade tip mach number limit and the advance ratio 
(uw) limit. The advancing blade tip mach number limit is set 
at M<0.8 to avoid compressibility effects on the advancing 
blade in forward flight. The advance ratio is set at p<0.4 
to avoid retreating blade stall at maximum forward speeds. 
One additional limit was incorporated to show the lower 
Reynolds number limit of 1.5 million. Fig. 7 depicts the 


constraint diagram developed using the above parameters. 
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Tip Speed Constraint Diagram 
Hummingbird 
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Figure 7 Main Rotor Blade Tip Velocity Constraint Diagram 


B. PRELIMINARY ESTIMATES 
1. Gross Weight and Blade Radius 
Ref. 3 outlines basic steps to quickly estimate the 
gross weight and rotor radius. Using these steps, the 


following calculations were made: 
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Fuel Required 
(0.5 lb/hp-hr) Piston Engine 
Mission time 60 minutes 
Fuel required = (1.0) (25) (0.5)= 12.5 lbs Gf Sie 
*This i8 approximately two gallons. 
Usable Load (UL) 
UL = crew + payload + fuel 
UL = 0 + 30 + 12.5 = 42.5 Ibs 
Assuming 1970 technology, Fig. 8 shows a useful load 
per gross weight factor of 0.4 [Ref. 3:p. 641]. 


Gross Weight = 90.6 lbs. 


Gross Weigni 


Q 
19S$0 1955 1960 1965 1970 1975 1980 1985 1990 
Year 





Figure 8 Historic Trends of Ratio of Useful Load to Gross 
Weight 


This number for gross weight falls within the range 
developed by the trends analysis, though it is on the low 


side for a payload of 30 pounds. Realizing that the design 
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weight of any design will grow, it is assumed that the 
higher weight is more reasonable; therefore, a 130 1b gross 
weight will be used for the following estimates. Based on 
this weight and the DL range from 1.0 to 2.0, the rotor 
radius ranges from 4.55 to 6.43 feet. This rotor radius 
range is validated by calculation of required blade area ag 
described in Ref. 2. Assuming a maximum forward velocity 


(V..) of 70 knots, Fig. 9 shows that blade area required is 


max 
5.5 square feet, which is 1.83 ft* per blade for a three 
blade system and 1.375 ft* for a four blade system. This 
Beansliates into a radius of 4.9 ££, a DL of 1.74 lb/ft* with 
a chord of 0.375 ft (4.45 in) for the three blade system, 
and a radius of 4.7 Bie a Devof ft. 8oeand a chord of 02292 Et 


(3.5 inches) for the four blade system. 
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Figure 9 Determination of Blade Area for New Rotor Design 
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Either of the above systems would be feasible for this 
design, but the four blade system does have two distinct 
advantages. First, the blade vibrations ina four blade 
system are less than those ina three blade system. Second, 
the aspect ratio (AR) of the three blade system is 13.1 and 
the AR for the four blade system is 16.1. The normal range 
for AR is 15 to 20. Any blade below that will have a lower 
blade efficiency due to tip losses and any blade above that 
could pose structural problems. 
2. One Hour Estimate 

In this estimate, there are two constraints: hover 
out of ground effect (HOGE) and high forward speed flight. 
The following is a list of the tentative performance 


requirements: 


Payload 20-30 lbs 
Crew O lbs 
Max Speed @ S.L. Zo knots 


Cruise Speed (0.9*V 


Max 


)}63 knots 

Vertical Rate of Climb 450 fom @ 4000 ft 95 deg. F 
Engine: One - Max continuous rating 25 BHP. 

The following calculation as detailed in Ref. 4 will be done 
using the previous gross weight estimation of 130 lbs. Fig. 
10 shows that for a weight of 130 lbs, the equivalent flat 
plate area (f,) would be approximately 0.8 £t*. [Ref. 5:p. 
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Figure 10 Parasite Drag Area vs Design Gross Weight 


a. Maximum Forward Velocity (V__) 


Assuming seventy percent of installed power is 
used to overcome parasite drag at high speeds, maximum sea 


level speed can be calculated. 


; : o 
Vin = 42*f (30 min fating eng) 13 (5) 


e 


Maximum BHP = 25 


Assuming Figure of Merit of 70% 


THP = 17.5 


V 114 knots 


Max 


Based on the assumption that the helicopter was to be 


derated to a power of 18 HP, the calculations would be 


changed as follows: 


9 


Maximum BHP = 18 


Assuming Figure of Merit of 70% 


TH? eel 
Vi. = 41* 112.6) eor ly 2) 
1! Se Oe oo) or 


max 


b. Rotor Sizing 


The preliminary size of the rotor is determined 
using the rate of climb (ROC) stated earlier. Additional 
power required to climb 450 fpm is approximately 10% of the 
power required for HOGE. 

ROC 450 fpm @ 4000 ft 95 deg. F (95% of 17 °33n 

95% of 17.5 tape = 16763 ue 

Pl = 1310) lby ferrocene aoe 

From hig. a1 the? Di~ =. 3 sri Rem ee aero le 

DL = Weight/(pi*R*2) = 3.5 

R = 3.44 ft 
Again assuming a derated engine, 

Seen Ore Aas nyo) S LZ wi) ays. 

Pi =) 3061/12 70min ee 

From Fig. 11 

Di = La75- 


DL 


Weight/ (pi*R*2) = 1.75 


R= 4.86 ft 
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Figure 11 Rotor Performance for Design Conditions 


c. Starting Point 


jhe Wreceding calicuiiabiegs provide a good point 


at which to begin more in-depth design. 


starting point is 
Gross Weight 
Empty Weight 
Payload 
Fuel Weight 
Rotor Radius 


Disk Loading 


In summation, the 
130 libs 
87.5 Ibs 
30.0 lbs 
ie Sieh Ds 
(R) = 4.86 


Oijs=-1.75 


ik 


3. Design from Scratch or Procure 

Based on this data, there were two choices to make. 
First, design a helicopter from scratch which would meet the 
needs of the Aeronautics and Astronautics Department. While 
this is definitely an exceptional learning process, the end 
result would leave the department with a paper helicopter 
and at least a year’s worth of manufacturing. The other 
path of choice was to look for a helicopter on the open 
market that closely approached the stated needs. The second 
alternative, being the most productive and time smart, was 
chosen. Though the possibility of locating a helicopter 
that would meet every need was slim, it was the best choice, 
realizing that the redesign of an existing, flying 


helicopter was much more frugal financially and time-wise. 
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Vv. HUMMINGBIRD ACQUISITION 


A. CAPABILITIES 

A helicopter of the correct size and payload capability 
was located at the GMP model helicopter company. The owner 
and designer, Mr. John Gorham, the lead engineer with 
Lockheed Corporation on the L-1011, had built ten 165 pound 
RPVs for the U.S. Army. Their original purpose was to be 
used as Soviet Hind-D recognition devices. This vehicle was 
nearly ideal for meeting the needs of the department. There 
were some shortfalls, which will be mentioned in the next 
section, but nothing of a critical nature that could not be 
redesigned or changed. The department purchased one of Mr. 
Gorham’s helicopters, with the desire to purchase a second. 
The initial RPH was named Hummingbird I. The purpose of 
buying two helicopters is to allow for concurrent HHC, 
NOTAR, and ALTOS research. One helicopter would be used in 
its original configuration for NOTAR and ALTOS research. 
The other would be converted into a platform capable of HHC 
research. Fig. 12 contains two pictures of the Hummingbird. 
Table I contains a list of the Hummingbird characteristics 


and capabilities as it was received from Mr. Gorham. 
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Table I HUMMINGBIRD CHARACTERISTICS AND CAPABILITIES 
Sm _ 


Hummingbird: 


Characteristics 


Weights 

Max Gross Weight 
Empty Weight 
Fuel capacity 


165 lbs 
15 TDs 
6.5 lbs 


Rotor Parameters Main Tail 


| Radius (R) i a blest 
Cintengen em, 6" 2.625" 
Solidity (sigma) 0.0061 Ae 

No. of blades (b) 2 3 

Tip speed 303 fps 241 fps 
Twist 25 0 

Hinge offset ration (e/R) }} 0.127 0.24 

Por tow 1 NACA 0012 NACA 0012 


Engines 


iy DS 

Number 

Maximum Usable Power 
Maximum Torque 


Westlake 342 Series 2100p, 


- 
25 BHP @ 7000 rpm 
25 ftlb @ 4000 rpm 





B. SHORTFALLS 

The Hummingbird satisfied many of the department’s 
needs; specifically, it was large enough to be developed 
into a quarter scale model, it had a 50 plus pound payload, 
and in addition, it had already proven itself in flight 
test. some shortcomings that fed to be 


There were, however, 


addressed. These shortfalls included the need for more than 
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a two blade hub, an RPV that possessed autorotational 
capabilities and could produce blade Reynolds numbers which 
would be comparable to full scale helicopters. 
1. Two Blade Rotor 

In order to have the capability to do HHC research 
work on the RPV, it was required that the helicopter have at 
least three rotor blades or more. The reason for this 
requisite is that the vibratory forces at the rotor blade 
root are produced by the (n-1)/rev, n/rev, and (n+1)/rev 
vibrations. The n/rev vertical forces and moments are 
transmitted to the fixed system at a frequency of n/rev. 
The (n-1)/rev and (n+1)/rev flapwise blade root shears 
result in n/rev hub pitching and rolling moments in the 
airframe. The n/rev flapwise blade root shears feed into 
the airframe as n/rev vertical forces. The (n-1)/rev and 
(n+1)/rev chordwise root shears produce n/rev airframe hub 
forces in the fore and aft and lateral directions. The 
n/rev chordwise root shears result in n/rev hub yawing 
moments. Choosing a rotor system with three or more blades 
therefore will not interfere with the primary 1/rev control 
inputs. [Ref. 6:p 19] This change alone created a 
Significant amount of work because it includes changing the 


main rotor hub and redesigning the rotor blades. 
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2. Lack of Autorotational Capability 
The original design never incorporated the ability 

for the Hummingbird to autorotate, and this was considered 
unacceptable for the NPS flight research vehicle. The 
bearings in the drive system do not allow the main rotor to 
freewheel should the RPV have an engine failure. Presently, 
should the engine fail it will cause a rapid decay in rotor 
rpm which will transform the RPV from a flying machine to a 
projectile. This is a critical redesign requirement which 
must be accomplished in order to protect the Navy’s 
financial and research investment in the Hummingbird. 

3. Low Reynolds Number 

In the HHC research, it is not enough to prove that 

the concept will reduce vibration; that result is 
sufficiently proved. The desire is to have the ability to 
test whether HHC also provides performance improvements. 
This comparison can only be accomplished on helicopters with 
similar Reynolds numbers. 

Reynolds number is the ratio of the inertia to the 
viscous forces on a volume of fluid (Eqn. 6). 


Re = pVve (6) 
I+ 
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The importance of Reynolds number in the comparison of 
geometrically similar bodies in incompressible flows can be 
shown using the Navier-Stokes equations. Eqn. 7 shows the 


dimensionless Navier-Stokes equation. 


Du’, 9p’, 1 g2,/, Dv__ ob, 1qi25) Dwi lee ae 
Dt! dz!) Re ae 


Dt! ax! Re —' pt! ay! Re 
These equation show that, given geometrically similar 
bodies, the equations of motion are identical for the same 
Reynolds number. The similar bodies includes surface 
roughness as well as shape. [Ref. 7:p. 304] 

The Hummingbird’s Reynolds number was on the order of 
0.9 million, and it needed to be in the range of 1.5 to 1.8 
meinen to be comparable to that of the OH-6A’S Reyvnolase 


number equal to 2.4 million. 
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Vi. BLADE DESIGN 


A. SPECIFICS 
1. Airfoil 

The airfoils that were available with the 
Hummingbird were NACA 0012 and NACA 0013. This type of 
airfoil is widely used in industry for numerous reasons. 
First, the airfoil is symmetrical and therefore there is no 
nose-down pitching moment which is associated with cambered 
airfoils. Second, it is a relatively thick airfoil which 
will provide an acceptable maximum lift coefficient. 
Finally, there is a vast amount of data available on this 
airfoil, which provides for much easier analysis. This is 
due to the fact that the NACA 0012 airfoil was selected by 
Mr. Sikorsky for the VS-300, the world’s first truly 
successful helicopter. It was also the airfoil of choice 
Hotealmostsall early helicopters, including more recent 
aircraft still in service, such as the Navy H-3. 

For improved performance, a more recent advanced 
technology airfoil worthy of consideration is the NACA 
23012. The characteristic drooped nose is an effective 
method of increasing the maximum lift coefficient. Also, at 
high lift it has lower drag than a similar six-series 


airfoil at low mach numbers. [Ref. 3:p. 388] 
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2. Twist 

The twist of a rotor blade enhances two main areas 
of the helicopter’s performance--hover performance and 
retreating blade stall--but it also produces increased 
vibration in forward flight. Its most notable adverse 
effect is to increase blade vibratory stresses and in this 
way reduce blade fatigue life. Built-in blade twist affects 
the radial variation of inflow angle from blade tip to blade 


root. Ideal twist is represented by Eqn. 8. 


(8) 


Hovering performance is enhanced with the addition of 
negative twist by creating a more uniform inflow 
distribution along the blade span. The larger the amount of 
twist, the closer it approximates an ideal twist 
distribution. Generally accepted values of main rotor blade 
twist are -8 to -14 degrees. Ideal twist yields the minimum 
induced loss for a given thrust. Retreating blade stall is 
delayed when twist is employed by unloading the tips, which 
reduces the tip angle of attack [Ref. 8:p. 57]. 

These regions of enhanced performance are not 
Critical areas for the Hummingbird. The Hummingbird will 
Operate at relatively slow forward speeds; therefore, 
retreating blade stall is not a concern. It would be 
attractive to have highly efficient blades in the hover 


regime, but with the present payload capability of 50 
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pounds, it is not critical. If these were critical areas 
for the Hummingbird, it might be worth the effort to 
manufacture new blades with twist; but since they are not, 
simplicity rules, and it was decided not to incorporate 
twist in the present rotor blade design. 

3. Tip Velocity 


Determining the tip speed (V,._) for the Hummingbird 


tip 
was a difficult matter because its effects were coupled with 
so many other areas of the helicopter’s performance. The 


constraint diagram from Chapter II is shown again in Fig 13. 


Bae design point 1s shown to be in the proximity of V_ of 
40 knots and es of 450 fps. V_,, was obtained from the 
Hummingbird original design. The means by which Viip Was 
determined will follow. Low tip speeds have the advantage 
of low noise and good hovering performance. "High tip 
speeds have the advantage of low rotor and drive system 
weights and high stored energy for autorotative entries and 
flares." [Ref. 3] In the case of the Hummingbird, Viip Was 
the variable used for a tradeoff study. Von had to be large 
enough to produce Reynolds numbers which could be compared 


with full scale helicopters, but small enough to provide a 


reasonable Figure of Merit (FM). 


a) 


Tip Speed Constraint Diagram 
Hummingbird 
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Figure 13 Constrain Diagram 


Bi. TRADEOFF 

The need to achieve a Reynolds number in the range 
mentioned earlier was the source of some consternation. If 
Reynolds numbers were not considered, it was quite easy to 
design a rotor blade that would fly efficiently and meet the 
department needs. The complication begins when the V ie 
Reynolds number, Figure of Merit and chord length 


requirements are all met at the same time. Reynolds number 


a2 


is a function of V,, and chord length as shown by Egqn. 9. 


up 
Re = 6400V,,,c (9) 


Fig. 14 and 15 show this relationship over a certain range 


Tradeoff Study 


Vtip vs Reynolds Number 
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Figure 14 Reynolds Number vs Blade Tip Speed 


of V,, and chord length. 
Figure of Merit is shown by the relationship in Eqn. 10. 
Induced power is the power required to overcome the drag due 


to the generated lift. Total power also includes the power 
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FM = ; (10) 





required to move the rotor blades through the air (profile 
power) and the power required to drag the fuselage through 
the air (parasite power). Fig. 16 and 17 show this 
relationship over the same Veo and chord length range. 

To link the Figure of Merit and Reynolds number together, 
lines of constant Reynolds numbers were plotted on the 
Figure of Merit versus chord length chart. Fig. 18 allowed 
for the tradeoff to be visualized simultaneously on one 
graph. The Hummingbird came with three sets of blades. One 
set was of radius 4.52 feet with a chord of 6 inches. The 
other two sets were of radius 5.0 feet with a chord of 6.5 
inches. To expedite the design process, the available 
blades were selected versus designing and manufacturing new 
ones. Since there was a requirement for at least a three 
bladed rotor, the second set of blades was chosen. 

The advantages of this choice are seen in the tradeoff 
study. Choosing a longer chord length allowed for a slower 
tip speed for the same Reynolds number. Since profile power 
is proportional to eee the slower tip speed increases the 


FM from 0248.00 70254. 
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Tradeoff Study 


Vtip vs Figure of Merit 
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Figure 15 Figure of Merit vs Blade Tip Speed 


C. STABILITY 

The pitch and roll damping is produced by the tilt of 
the tip path plane, which lags behind the motion of the 
shaft by an amount that is proportional to the rate of pitch 
or roll. Therefore, the aerodynamics on the blade causes 
the tip path plane to tend to stabilize itself in an 


equilibrium position with respect to the shaft. "If the 
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aerodynamic and inertia flapping moment are equated, the 
following results for the angular displacement of the rotor 
plane with respect to the shaft per unit tilting velocity of 


the shaft is obtained for the hovering case:"[{Ref. 8:p. 275] 


Tradeoff Study 
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Figure 16 Reynolds Number vs Chord Length 
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The quantity 16/(y2) can be understood more easily by 


examining its physical interpretation as follows: The thrust 


a6 


vector lags the rotor shaft by a time constant of 16/(y¥N) 
seconds if the rotor shaft is tilted with any constant 
angular velocity. Therefore, the larger the time constant 


(16/(y)), the greater the system is damped. Lock number 


Tradeoff Study 
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Figure 17 Figure of Merit vs Chord Length 


(y) relates the inertia and aerodynamic characteristics of a 


blade and is shown in Eqn. 12. 


y _ cpar* (12) 
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Figure 18 Figure of Merit and Reynolds Number Tradeoff 


The stability of the Hummingbird is a critical area that 
will require additional research. The two proposed 
modifications to the helicopter significantly alter the 
stability of the RPV. Eqn. 11 shows that rotor speed and 
Lock number will affect the damping of the rotor. The rotor 
blade is not changed; therefore, the Lock number will remain 
the same. The rotor speed will be increased from 550 rpm to 
759 rpm, thus decreasing the effective rotor damping. Ref. 


8 states that 
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In addition to the effects of rotor speed, rotor damping 
may be increased by the use of devices that act upon the 
control system in such a manner as to increase the 
displacement of the rotor from its trim position due to 
a given rate of roll or pitch. An example of such a 
device is a rate gyro that would apply opposite control 
by an amount proportional to the rolling or pitching 
velocity of the helicopter. 


The Bell-Hiller stabilizer bar hub configuration is such a 
device. Therefore, changing from that type of hub to an 
articulated hub will significantly decrease the rotor roll 
and pitch damping and thus the helicopter’s handling 
characteristics. For the 3-bladed Hummingbird, rate gyros 
may be required in order to have satisfactory flying 


qualities. 


D. BLADE DYNAMICS 

1. Introduction 

tt isevery important to consider the vibratory 

effect on the main rotor blades. High vibrations result in 
high vibratory shear at the rotor hub, which results in high 
fuselage vibrations. The high vibrations also cause 
vibratory shear stress, which reduces the effective life of 
a rotor blade. The vibratory resonance response also 


affects dynamic stability of the rotor system. 
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2. Main Rotor Blades 

The following analysis was done on two different 
rotor blades. The first was termed the "heavy" blade. It 
had a length of 4 ft 6.25 inches, chord of 6 inches, 
thickness ratio of 12 percent and a weight of 5.18 lbs. The 
second blade, termed the "light" blade, has a length of 5 
ft, chord of 6.5 inches, thickness ratio of 13 percent anda 
weight of 4.66 lbs. The length given for the blades above 
does not include the 8 inch blade offset of the rotor hub, 
but it was included in the rotor radius when the analysis 
was performed. 

In order to analyze each blade, they were broken up 
into segments: the heavy blade into 17 segments and the 
light blade into 24. Fig. 19 and 20 show a schematic for 
the two blades and how the blades were divided into 
segments. The analysis for each blade was conducted with 
the exact same technique, and therefore the discussion will 
be limited to the heavy blade. The following quantities 
were required to perform the blade analysis: section radii, 
chord, segment width, segment volume, segment weight, and 
area moment of inertia (Ixx), and material modulus of 
elasticity. The section radii was the distance from the 
center of rotation to the center point of the segment. The 
chord was constant for the blade until the small taper at 


the root. The segment width was constant with the exception 
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of the tip segment, which was allotted any leftover blade. 
Fig. 21 shows a cross section of the blade, which was used 
for area and moment of inertia calculation. For 
calculations, the airfoil cross section was broken up into 
three sections and approximated as follows: (1) the nose 
section by a parabola; (2) the center section by a box; and 
(3) the tail section by an isosceles triangle. 

The cross sectional area of each section was 


calculated using standard geometric formulas in Eqn. 13. 


Parabola: A = 420 


Box: A= bh Triangle: A= Sbh aS 


The area moments of inertia for the box and triangle were 


Meavy Blade 
Wuight = 5.18 lbs 





Figure 19 Heavy Blade 
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Light Blade 


Weight = 4 66 lbs 
Co a | 
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Figure 20 Light Blade 


calculated using Known formulas, while integration was used 


to determine the Ixx of the parabola. (Eqn. 14.) 
Box: Ixx =—- bh? Triangle: Ixx = bh} Parabola: Ixx | y*da 
de, IP A 


(14) 
The parabola was approximated by the equation x = 6.67y’. 
Using the calculated cross sectional area and the 
segment width, the volume was calculated. The segment 
volume was then used to determine the segment weight. The 
assumption was made that the rotor blade was made of a 
homogeneous material - spruce. The only variation in these 


calculations was where the blade was tapered very near the 
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Hummingbird Airfoil Section 





Figure 21 Main Rotor Blade Cross Section 


root. This section was no longer an airfoil section, but an 
inch thick tapered block. The segment weight and Ixx of the 
root section also included contributions from the blade grip 


assembly. Table II shows the data for each blade. 
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Table II BLADE STATION DATA 
ar 


Heavy Biade Radius = 452 ft 
Weight 5 18 Ibs Chord = 6i1n 
Uc = 12% 


Section Radil Chord Xsec Ar Seg Wd Volume Seg Wt. ixx Delta Cf 
(in) (in) = (in* 2) (in) = (in * 3) (lb) {in* 4) (ftib) 
80 


1 59.1 168 6 25 1050 0.459 0 097 230 65 
2 545 60 168 300 5 04 0 220 0 097 49 03 
3 515 60 1 68 3 00 5 04 0 220 0097 46 33 
4 485 60 168 3 00 5 04 0 220 0 097 43 63 
5 455 60 1 68 3 00 5 04 0 220 0.097 40 93 
6 42.5 60 168 300 5 04 0 220 0 097 38.23 
ri 39 5 60 168 3 00 5 04 0 220 0.097 35 53 
8 365 60 168 300 5 04 0 220 0 097 32 84 
oi 335 60 168 3 00 5 04 0 220 0 097 30.14 
10 30 5 60 168 3 00 5 04 0 220 0.097 27.44 
1 275 60 168 3 00 5 04 0 220 0 097 24.74 
12 245 60 1 68 3 00 5 04 0 220 0097 22 04 
13 215 : 60 1 68 3 00 5 04 0 220 0 097 19 34 
14 185 60 1.68 3 00 5 04 0.220 0 097 16 64 
15 155 5.4 5 40 3.00 16 20 0 708 0.450 44 82 
16 125 50 5 00 3 00 15 00 0 656 0.417 33 47 
17 95 38 375 3 00 11.25 0742 0 808 28.77 
Light Blade Radius = 50ft 
Weight = 466 ibs Chord =65in 
Uc = 13% 


Section Radli Chord Xsec Ar Seg Wd Volume Seg Wt. xx  Deita Ct 
(10) (in) = (in * 2) (in) = (in * 3) (ib) fin* 4) (fib) 


1 65 6 65 3 37 48 16 18 0353 0 478 151 24 
2 62 0 65 337 24 8 09 0177 0.478 35.74 
3 59 6 65 3 37 24 8 09 0.177 0.478 34.35 
4 $72 65 337 24 8 09 0.177 0.478 32.97 
5) 548 65 337 24 8 09 0.177 0.478 31.59 
6 92.4 65 337 2.4 8 09 0.177 0.478 30.20 
7 50 0 65 337 2.4 8 09 0177 0.478 28 82 
8 476 65 337 24 8 09 0.177 0.478 27.44 
9 45 2 65 3.37 24 8.09 0177 0.478 26.05 
10 428 65 3 37 24 8 09 0.177 0.478 2467 
11 404 65 337 24 8 09 0177 0478 23.29 
12 38 0 65 337 2.4 8 09 0177 0.478 21.90 
13 356 . 65 337 2.4 8 09 0.177 0.478 20.52 
14 33 2 65 337 2.4 8 09 0177 0.478 19.14 
15 30 8 65 337 24 8 09 0.177 0.478 17.75 
16 28 4 65 337 2.4 8 09 0.177 0.478 16.37 
7 260 65 3.37 24 8 09 0.177 0.478 1499 
18 236 65 337 24 809 0177 0.478 13 60 
19 212 65 337 24 8 09 0.177 0.478 12.22 
20 188 65 337 2.4 8 09 0177 0478 10 84 
21 164 65 337 24 809 0177 0.478 945 
22 140 38 5 76 24 13 82 0 302 0.480 13.79 
23 116 48 4 80 24 1152 0 252 0 400 9.52 
24 92 42 420 24 1008 0 470 0845 14.1] 
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3. Myklestad Determinant Method 
There are many methods for analyzing vibrations ina 

rotor blade. Some of these are 

® Rayleigh-Ritz 

® Holzer 

® Myklestad-Prohl 

® Matrix iteration 
The Myklestad method, which determines the natural 
frequencies and modes of the blades, was the method of 
choice. It is especially well suited for vibration analysis 
of rotor and turbine blades. Dr. Nils Myklestad was a 
consultant to Bell Helicopter Company for many years. The 
proper analysis requires considering the coupled flapwise- 
edgewise-torsional response of the blades. Because the 
edgewise and torsional stiffness is much greater than the 
fFlapwise stiffness and thus has much high frequencies, only 
the flapwise responses of the blades were considered. 

The uncoupled flapwise Myklestad system is shown in 

omen 22. 
From this diagram the equations of equilibrium for this n‘ 


element can be written. 
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Figure 22 Flapping Blade Diagram 


Centrifugal Force: 


= 2 
Ty = Ty +m 2% ry 


Shear: 
Sy=Syiy tM y @ : Yn 


Moment: 


My = Mya tSyoy er rey (Yui -Yw) 
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(15) 


(16) 


(17) 


Slope: 


Oy = Oye (14T Vasnst) oy, Tovat_g  Unsne (18) 
N N +1 N +1 QEL +] EL N+1 QEIL 


Deflection: 


3 2 3 
Tye Oyi 1 n,n+1 Mal n,ast_ Syart n,n+i (19) 


= =O < aa 
Yn Yw+t N n,n+i 3ET Pie 3EIL 


The natural frequencies of the system are obtained 
by assuming a frequency (Q) and then calculating the 
centrifugal forces, shears, moments, slopes, and deflections 
of each blade segment from the blade tip to the blade root. 
The boundary Pere at the blade tip are 
| Ss, =M = 0; 6, = 86,;; y, = Y,; 
pememe TOOL Of the blade, the two unknowns 9, and y, are 
carried along and the equilibrium equations are written as 
Follows: 

So = Asy; + B,O, 
o = Avy; + By8, 
o = Agy; + BO, 

Y, = AY; + Bye, 


In the preceding equation, the coefficients A and B will be 


S 
T 


cep) 
ll 


calculated as the process proceeds down the blade. They are 
functions of mass and stiffness properties of the blade, the 
rotational speed, and the assumed frequency. The boundary 


conditions at the root of the hinged blade are 
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Or 


0 


AyY; + By®, 
AVY5) BO p = a0 
These equations can be written in matrix form. The 
frequencies that satisfy the equation (when the determinant 
equals zero) are the natural frequencies and can readily be 
determined. 

A Fortran code, Appendix A, was written to solve the 
Myklestad equations and determine the flapwise natural 
frequencies of both Hummingbird rotor blades. Fig. 
23,24,25, and 26 show a graphical output of the Myklestad 
determinant for both blades. 

The natural frequencies, where the line crosses the 
x-axis, are clearly seen. Appendix B contains a numerical 
output of the same data for the different natural 
frequencies. 

A means to verify the accuracy of the code’s ability 
to determine the natural frequencies was to check the ratio 
of the determinant coefficients. A, and A, are the 
displacement and slope Peele ea eeci. ively for the 
rotor blade root section in response to starting boundary 
conditions at the blade tip of slope equal zero and 
deflection equal one. B, and B, are similar to iy and A, 
except the starting boundary conditions are slope equals one 


and displacement equals zero. Knowing that the determinant 
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Figure 23 Myklestad Determinant (Heavy Blade) 
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Figure 24 Myklestad Determinant (Heavy Blade) 
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Figure 25 Myklestad Determinant (Light Blade) 
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Figure 26 Myklestad Determinant (Light Blade) 


equals zero at the natural frequency, the ratio of Ny) =. 
should equal -A,/B,. Using data in Appendix B for the heavy 


blade’s first natural frequency, the accuracy is seen. 


ist Natural freq. = 107.1522 rad/s 
A= 7.48560 BS -469.659 
A = 1.79121E+04 B= -1.12384E+06 


2X2, = 1595837 h-02 Ae/Be a 1 SoS 550E 02 


S2 


It is interesting to note that the natural 
frequencies are much higher than those of conventional 
helicopters. This is attributable to the fact that model 
rotor blades are much shorter and made of solid spruce, 
which makes them considerably stiffer than full scale 
blades. 

4. Blade Modes 

The blade modes were used to verify that the 
calculated frequencies were in actuality the ascending 
natural frequencies. Appendix C contains a Fortran code 
written by Lt. M. Avila that was modified for use on the 
Hummingbird rotor blades to obtain the rigid and bending 
flapwise mode shapes. Fig. 27 shows the mode shapes for the 
heavy blade while Fig. 28 shows the mode shapes for the 


light blade. 
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Figure 27 Heavy Blade Mode Shape 


5. Southwell Plot 
The Southwell plot, also called a Fan plot, isa 
resonance diagram of the individual blade. It is used to 
depict the flapwise modes in relation to the integer 
multiples of the rotor rpm, commonly referred to as 1P, 2P, 
3P etc. For a good design, effort should be directed toward 
keeping the blade natural frequencies away from the Ws) OE 


etc. forcing frequencies, since these are the primary 
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Figure 28 Light Blade Mode Shape 


excitations. The Fortran code used for the Myklestad 
determinant was run for rotor speeds ranging from 0 to 90 
rad/s. The results of these runs are shown in Table III. 
The area of concern is within plus or minus 6 rpm of the 
operating rpm. Fig. 29 and 30 show the Southwell plot for 
the heavy blade and light blade respectively. 

The rigid mode always lies very near the 1P line. 
The flapwise bending modes for both blades are much higher 


than that for full scale blades. Again, this is attributed 


5) 1S) 


Table III BLADE FLAPWISE RESPONSE 





Heavy Blade 


Bend Bend Bend 
Pf tepm | cpm | cep | (opm) | tepm 
fo [1.29 {1758.9 | 5757.3 [1.19 | 2876.2 | 
2093.8 


2945.7 
30 3030.0 
30 3244.2 
50 3285.4 
3449.2 


Wis 1023 a 3.222 JZLEC9 gy oe. 7 3824-3 


9 
90 |eso.4 | 1159.3 | 3535.2 | 7608.9 








to the greatly increased stiffness. The second and third 
modes for the light blade were higher than those for the 
heavy blade. This is accounted for by the fact that 
although the light blade is larger, the lower mass more than 
offsets the length by increasing the natural frequency. The 
heavy blade Southwell plot show that there will be no 
vibratory resonance problems with any of the modes. The 
first flapwise bending mode of the light blade lies very 
near the 5P line at the operating frequency. The 5P lies at 
62.56 Hz while the second natural frequency lies at 63.74 


Hz. The close proximity of the two frequencies should not 
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be a problem for the Hummingbird, as the critical 
Frequencies occur at the 2P, 3P, and 4P for a three blade 
rotor system. For a full-scale equivalent 3-bladed 
helicopter, this near 5P resonance might need to be 
monitored due to possible high amplification of 5P flapwise 
blade stresses. In the case of the Hummingbird, the 
relatively short blades are sufficiently overdesigned to 
avoid this problem. 

Discussion of blade resonance results for the RPH blades 
would not be complete without comparing these values with 
values obtained for full-scale helicopters. Here for the 
light and heavy blades, we find the first rigid mode at 
1.262 and 1.35Q2 respectively. For a uniform articulated 
blade this would be considerably closer to 1P at 1.029.. The 
first flapwise bending mode for an articulated fullscale 
blade would be in the vicinity of 2.52, whereas analysis 
shows in this case that we have much higher frequencies of 
5.0392 and 4.2592 for the light and heavy blades respectively. 
Similarly, the second flapwise bending mode for a full-scale 
articulated blade is about 4.792 to 4.992, whereas for the 
heavy blade, the second flapwise mode occurred much higher, 
at 9.540. This is even higher than where we would normally 
expect the frequency of the third flapwise mode to occur for 
a full-scale rotor. Normally the third mode resonance 
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Figure 29 Southwell Plot (Heavy Blade) 
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Figure 30 Southwell Plot (Light Blade) 
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VIZ PERFORMANCE ANALYSIS 


A. POWER CALCULATIONS 

1. Power Required 

Using the parameters established in the blade design 

section and the procedures outlined in Ref. 9, power 
estimations were calculated. There exist two regimes of 
operation that establish the power required limits. The 
first is hover out of ground effect; the second is high 
speed flight. Since the Hummingbird would not operate in 
the high speed flight regime, the HOGE was established at 
the limiting flight condition. Two hover conditions were 
evaluated: sea level standard day and 4000 ft MSL 95 degrees 
F. As alluded to earlier, there are three basic power 
calculations that must be accomplished. They are induced 
power (P.), the power required to produce lift; protile 
power (P,), the power require to push the blades through the 
air; and parasite power (P,), the power required to move the 
rest of the fuselage through the air. The following 


formulas were used to calculate these quantities. 


ak 
1.5 2c.) 2 
p, = if HP: p-1- 22 (20) 
V2pA B 550 b 
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Cy, POA Vigoss 


= 4400 
i 3 
—pf.V; 
pee SP ze 
“ 5 518. 
fo— inrust B— lip boss Factor b = Number of blades 
A = Rotor Disk Area V, = Forward Velocity 
mes tQuivealent: Eldt plate area 


Using the above formulas, the total power required to hover 
out of ground effect was calculated. Induced power was 4.95 
HP, while profile power was 4.4 HP. The parasite power is 
equal to zero in the hover case. To calculate the power for 
the 4000 ft 95 degree F case requires a density correction. 
The power required to hover in ground effect (HIGE) was also 
calculated. The induced power ratio equation, based on the 
a ratio of the rotor height above the ground over the rotor 


diameter, is used for this calculation, and is shown below: 


( IGE) 


P, 4 3 
pee = 1276 (2) +0.7080 (2) -1.4569 4 
P, (OGE) D D 


tb lb 


2 h 
) 41434232 (—) 40.51 
D 
(23) 


The rotor height above the ground (h) equals 2.5 feet. 
Using the light blades for these calculations, the ratio of 
h/D equals 0.22. It can be shown that the induced power 
required to hover IGE is 3.7 HP. These calculations are 
solely for the main rotor system. The power required for 


the tailrotor will be included in later calculations. The 


Ge 


final flight regime that must be considered is forward 
flight. The profile power in forward flight is shown by 


Eqn. 24. 


Eorlight i. ever 1 t4 ssp) (24) 


This equation does not include high speed effects such as 
compressibility or retreating blade stall, but as stated 
previously, these will not be factors for the Hummingbird. 
Power required calculations for velocities ranging from 


hover to 100 knots are shown in Fig. 31. 
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Figure 31 Power Required Profile 
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Tailrotor power is included in these calculations, assuming 
that it is nine percent of main rotor in the — 
condition, decreasing to three percent at mid-range 
velocities. Appendix D contains the spreadsheet calculation 
for all power required calculations. 
2. Power Available 
The power available curve, Fig. 32, was constructed 


from data obtained from Ref. 10 and 11. An assumption of 


near linearity over the entire RPM range was made. 
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Figure 32 Power Available vs RPM 


ous 


Fig. 32 will be used again when determining the derated 
power that will be use in conjunction with resizing of the 


main rotor drive gear for obtaining desired rotor RPM. 


B. RETREATING BLADE STALL 

A phenomena that is particular to helicopters is 
retreating blade stall. It occurs at high forward speeds 
when the retreating blade is unable to produce the lift 
required to maintain equilibrium with the advancing blade. 
This phenomena is due to the great difference in local 
velocity over the airfoil. The advancing blade encounters 
the velocity of Viip Plus the forward velocity. The 
retreating blade encounters Veo minus forward velocity. At 
high forward speeds, the velocity over the retreating blade 
becomes very small, requiring a very large angle of attack 
to produce the needed lift for equilibrium. Also, the flow 
over the inboard part of the retreating blade is reversed, 
flowing from trailing edge to leading edge. As speed of the 
helicopter increases this reversed flow region moves out 
radially along the blade. Typically at cruise speed as much 
as 30 percent of the blade is experiencing reversed flow. 

The total angle of flow can be determined through either 
of two equations. The first equation is based on the 
average lift coefficient, which is lift curve slope (a) 
times the angle of attack, which will be expressed in terms 


of blade loading (C,/o). The second equation is based on 
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the flapping angle of the blade. Both of these equations 
aoe ice tons Or cOlUceeingemmaten angle, 9,, longitudinal 
aaelnempitch angle, U5, geometric angle of Ewist, 0, and 
inflow ratio, A. For each of these two equations, the four 
parameters are multiplied by coefficients which are 
functions of the tip loss factor (B), and the advance ratio 
(wu). The relationship for the first equation is given by 
Eqn. 25 with the accompanying coefficients given in Eqn. 26, 


Pweee2o and 29. 


Ze 


= iO nO. tO, 1, (25) 
Ea Oo BatO. ot) (26) 
ee 042 nes Oe Sita D) (27) 
1 eee A's) (sia see Uy (28) 
TOS eee a0 252°) (29) 


The second relationship is for the longitudinal flapping 
Seerttci1ent, a,, which is assumed to be zero or very close 
to it. The second set of relationships are given in Egn. 


Sop si, 32, 33, and 34. 


a, =0 = AA,,+0,4,,+074,,+0,A,, (30) 
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Ae 2 (31) 
**  B2(B2-0.5p2) 
8uB 
an pRB (32) 
‘2-3 (B2-0. 5p?) 
2 
A,; = __ 2pB" (33) 
BeOS 
2 2 
ad Bote oe (34) 
B*-0.5p? 


Solving the two equations simultaneously in terms of the 
eight coefficients, the collective pitch angle and the 
longitudinal cyclic pitch angle can be determined. The 


resulting equations are as follows: 


at aia (35) 
2 
2 Cr 
Aya th Ar Te-AreT1) +8 9(A,3 Ty -Ar 275) (36) 


0, = 
A121, Ay, T2 
Using these equations, approximations for the angle of 
attack of the advancing and retreating blades can be 
calculated. The equations for the angle of attack are given 


by Eqn. 37. 
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Goa 7 Og 7e510,t 





A 
Ty Goce aro omy (37) 


Appendix E contains the spreadsheet that was used to 
calculate the angle of attack for three different blade 
twist conditions, 0, -4, and -8 degrees of twist. The minus 
sign indicates that the twist angle decreases from the tip 
of the blade to the root. Fig. 33, 34, and 35 show the 


results of these calculations. 
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Figure 33 Retreating Blade Stall, 0 Twist 
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Figure 34 Retreating Blade Stall, -4 deg Twist 
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Figure 35 Retreating Blade Stall, -8 deg Twist 
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According to Ref.8 the blade reaches its operational limit 
when the angle of attack of the tip of the retreating blade 
exceeds the stall angle of attack (AOA) of the airfoil 
section by four degrees [Ref. 8:p. 258]. For the NACA 0012, 
this AOA would be 16 degrees. There are three areas of 
interest on these figures: airspeeds which require less than 
12 degrees AOA where no stall will be encountered, airspeeds 
which require between 12 degrees and 16 degrees AOA where a 
moderate amount of stall is present, and airspeeds which 
require 16 degrees or greater AOA where the severity of the 
stall prohibits flight. Two lines are annotated on the 
above graphs to show the velocity for this stall area. The 
limiting velocities for the three conditions are 64 & 74, 66 
& 76, and 68 & 77 knots respectively. The result agrees 
with theory, because the increase in twist should increase 


the stall envelope, which the figures demonstrate. 
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Vit T MODIFICATIONS 


A. LANDING GEAR 

The rear landing gear of the Hummingbird has a narrow 
tread width which is considered to be too narrow in view of 
takeoff and landing stability. It also has no shock 
absorption capability, which is critical for avoiding ground 
resonance. A study was conducted to determine the most 
advantageous route to take to overcome this challenge. The 
options were to leave it as it is, replace it with a skid 
configuration, or a strut configuration. The landing gear 
was redesigned with a strut configuration to overcome both 
of these problems. The tread width was widened from 17 
inches to 27 inches and an oleo was incorporated. The 
spring constant of the oleo was chosen so that the 
Hummingbird could survive a fall from ten feet. Ten feet is 
approximately one rotor diameter, the height to hover out of 


ground effect. 


B. DRIVE TRAIN MODIFICATIONS 
1. Engine to Main Rotor Gear Ratio 
Because of the requirement for the main rotor blades 
to obtain a prescribed Reynolds number range, the rotor rpm 
had to be increased to achieve the tip velocity required. 


Before the gearing ratio was chosen, the decision to derate 


Teak 


the engine had to be made. It was shown earlier that for 
the Hummingbird’s range of operation, there would always be 
excess power available. Based on this fact and the desire 
to extend the life of the engine, it is recommended that the 
engine be rated to operate at 18 HP vice 25 HP. Fig. 32 
shows that for this horsepower an engine rpm of 5000 is 
required. The present gear ratio is 10:1 and would have to 
be lowered to 6.6:1 to achieve a rotor rpm of 759. 
2. Freewheeling Unit 

The Hummingbird, as stated earlier, has no ability 
to autorotate. Because of this operating limitation, the 
Hummingbird is at great risk of loss to an inflight engine 
failure. It was decided that the most effective way to 
solve this safety hazard was to incorporate a one-way 
bearing in the gear that would be changed for the gear ratio 
change. After examination of the drive train, the last gear 
in the drive train was selected as the most effective point 
to make the modification. This point was selected for two 
reasons. First, because it was the last gear in the drive 
train, the rotor would have the least mechanical resistance 
during an autorotation. This point is critical because 
every rpm gives that much more cushion in the touch down 
phase of an autorotation, thus increasing the chances of a 
successful autorotation. Second, this gear has easy 


accessibility and would simplify the modification process. 
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Ix. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The goal of this research effort was to design a 
remotely piloted helicopter (RPH) with the capabilities to 
meet the expanding needs of the Department of Aeronautics 
and Astronautics helicopter flight test program. These 
needs have grown rapidly over the past two years. J. G. 
Scott stated in his thesis that the "requirements were: (1) 
a four-bladed rotor head; (2) payload capacity of 
approximately 15 pounds; and (3) a total system cost of no 
greater than $10,000." [{Ref. 12] Numerous research 
requirements have arisen since that time, which have 
expanded the envelope of RPH requirements. Lt. M. Borno’s 
quarter scale NOTAR research requires a platform of greater 
size and power than is available with any of the 
department’s current assets. Continuing HHC research 
requires an RPH that will produce blade Reynolds numbers in 
Biemprange of 1.5 to 1.8 million in order to validate HHC 
performance enhancements. Based on these two research 
requirements and the hindsight of how quickly these 
requirements change, it was imperative to develop an RPH 
that could meet the existing needs and would be flexible 


EPougnweo mect future needs. 
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Upon completion of the initial design phase it was 
determined that two paths were available: (1) complete the 
detailed design in house; or (2) search the open market to 
find an RPH capable of meeting the existing requirements. 
The second path was chosen and resulted in procurement of 
the Hummingbird. The Hummingbird possesses the capabilities 
to meet existing and future requirements with some 
modifications. The NOTAR requirements are fulfilled in the 
Hummingbird’s present state, but the HHC requirement 
dictates some modification. These modifications include a 
new three blade main rotor hub to be used with the existing 
rotor blades and a change in the engine to main rotor 
gearing ratio. Two other modifications that would be of 
great value from a safety standpoint are incorporating a 
freewheeling unit in the engine drive system and replacing 
the existing landing gear. 

The performance analysis shows that with the main rotor 
hub change, the Hummingbird will be able to achieve a 
Reynolds number within the range required for aerodynamic 
performance comparison with the OH-6A. This change will 
reduce the Figure of Merit, but due to low disk loading, 
will have little effect on the Hummingbirds power loading. 

The analysis of the Hummingbird also shows that it will 
have an excess power of nearly 10 HP for HOGE and therefore 
can be operated at a derated power. This would reduce wear 


on the engine and thus enhance its life. The decision to 
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derate the engine should be made prior to changing the 
engine to main rotor gearing, as it will effect the 
operating rpm of the engine. 

The design improvements; modified landing gear anda 
freewheeling unit, will provide the capability to safely and 
successfully recover the Hummingbird from nearly any 
malfunction over a large flight regime. 

The overall analysis of the Hummingbird is that it is an 
ideal platform for subsequent research. It has the 
capabilities to meet the existing needs and the flexibility 


to meet any foreseeable future needs. 


B. RECOMMENDATIONS 

There is almost unlimited potential for the 
Hummingbird’s use in further RPH research. Two items that 
need to be undertaken as soon a possible are the design and 
manufacture of a three blade hub and incorporating a 
freewheeling unit in the engine drive train. The first is 
necessary before any HHC work can be undertaken on the 
Hummingbird. Measurement could then be started on the 
control system freeplay and torsional constants. The ground 
work is laid out in Ref. 12 and 13. The second item is 
primarily for safety reasons and to help protect the 
investment of the department. 

The design of a new main rotor hub could be a thesis in 


itself. A great source of expertise that could be of help 
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is Mr. Art Phelps of the Aerostructures Directorate at NASA 
Langley. His directorate owns an RPH identical to the 
Hummingbird, and it would be of value to develop the bonds 
of cooperation that have been established with them. 

Two other areas where valuable research could be done is 
in the design of an improved main rotor blade and the 
incorporation of a new landing gear design. An improved 
main rotor blade should be designed to optimize airfoil 
cross section and blade twist. The manufacture of new 
landing gear should be a moderately easy job which would 
greatly enhance the survivability of the Hummingbird. 

The final recommendation is made in light of all the 
current and possible future research that would hinge on 
the Hummingbird. It is highly recommended that a second 
Hummingbird be purchased. The reason for this is that 
concurrent research on the Hummingbird will inevitably 
require alteration that would adversely affect one or both 
of the areas of research. A case in point is the hub 
redesign and how it will affect the NOTAR research. The 
NOTAR research demands that the rotor operating parameters 
be constant. This demand would not be feasible in the 
tradeoffs that occur in the design process. Therefore, 
having a second helicopter would allow for numerous research 
projects to be accomplished concurrently without adversely 


affecting one another. 
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APPENDIX A: MYKLESTAD FORTRAN CODE 


THIS PROGRAM CALCULATES MYKLESTAD DETERMINANT FOR 
DETERMINING THE NATURAL FREQUENCY OF A ROTOR BLADE 


THIS PROGRAM IS FOR AN ARTICULATED HUB ONLY 


QANAAN 


THE PROGRAM IS SET UP TO TAKE A BLADE OF 17 SEGMENTS. 


IMPLICIT DOUBLE PRECISION (A-H,L,M,0-Z) 

DOUBLE PRECISION T(1:17),SIFEAR(1:17) ,MOM(1:17),SLOPE(1:17), 
$¥(1:17),R(1:17) ,WT(1:17),DELTA1, XI (1:17) ,NR,NRSQ, F,WSQ, 
%$L(1:17),W,S1,S2,83,Y1,¥2,Y3,¥4,MOM1,MOM2,8S, BM, DELTA2, 
$BSLOPE, BY,AS,AM,ASLOPE, AY, XPLOT(1:17) , YPLOT(1:17) 


CHK RK KK KKK KK KKK KKK KKK KK KKK KEK EE KEKKKEEKEKEKKKKKKKKKKKKKKKKKKKEK 


C VARIABLE DEFINITION 

C 

eT - CENTRIFUGAL FORCE 

C SHEAR - SHEAR 

C MOM - MOMENT 

C SLOPE - SLOPE 

Cry - DEFLECTION 

“eee eee eee ee eT eee Cee eee SESE SESE SS SSC SCS SCCCSOCSCOCSSCS CICS TT 


c : 
C THIS OPENS THE DATA FILE FOR OUTPUT “WOODPTS.DAT" 
Cc é 


OPEN (UNIT=35,FILE='WOODPTS. DAT’ ,STATUS=’ UNKNOWN ) 
C 


CREKHRRKKKKEKKEKKERKEKEKKEKKEKRKEKRRKRKRRAKRKRKKKRKKEKKKEKKKKKKKKKKKKKKkKkkhKEK 


C 
C THESE ARE THE INITIAL PARAMETERS TAKEN FROM THE BLADE DATA 
C THEY ARE FROM THE BLADE TIP TO THE BLADE ROOR 
C 
C SEGMENT SPACING (IN) 

DATA L/6.25,16*3.0/ 
C SEGMENT RADIUS (IN) 


DATA R/59.1,54.5,51.5,48.5,45.5,42.5, 39.5,36.5,33.8,30.5, 

$27.5,24.5,21.5,18.5,15.5,12.5,9.5/ 
C SEGMENT WEIGHT (LBS) 

DATA WT/.459,13*.22,.708,.656,.742/ 
C SEGMENT FLAPWISE MOMENT OF INERTIA, Ixx (IN*4) 

DATA XI/14*.097,.45,.417,.808/ 
CHI KKK KKK IK KKK KK KEK KKK KKK KEK KEK KEKE KEKE KKEEKKKEKEKKEKKKKKEKKEKKEKER 
C ALLOWS FOR DIFFERENT NR INPUTS WITHOUT RECOMPILING 
Cc 

PRINT *, ‘INPUT ROTOR RPM OMEGA’ 

READ *, NR 

NRSQ=NR*NR 

E= 1305000 ! MODULUS OF ELASTICITY OF SPRUCE WOOD 


CR IK KK KKK KKK KK KI KKK KK KKK KKK KKK KKK KEKE EKKKEKKEEKKKKEK 
C THIS LOOP FOR F MAY NEED TO BE MODIFIED FOR THE LOWER FREQUENCY 
C RUNS AS YOU DO NOT OBTAIN ENOUGH DATA 
Cc 

DO 400 F=.1,30,.1 

W=(NR*F) 

WSQ=W*W 

I=1.0 


Cl 


CR III TOI RIOT TOTO TOI IOI TOTO TOR TOR TTT aK 
Cc 


C INITIAL CONDITIONS 
Cc 

SLOPE(1)=1 

Y(1)=0 
020 SHEAR(1)=0 

MOM(1)=0 
G 
CH Hee tHe KKH te HK THK Ki KKK KK HTK HK KK TK TT KK KKK KKH KKKKKKKKK 
é 


C MYKLESTAD METHOD 


T(1)=WT (1) *NRSQ*R(1) /(32.174*12) 
DO 30 J=2,17 
T(J) =T(J-1) +( (WT (J) *NRSQ*R(J)) /(32.174*12)) 
S1=SLOPE(J~1) *(1+((T(J-1) *L(J) *L(J)) /(2*E*XI(J)))) 
S2=MOM(J-1) *L(J) /(E*XI (J) ) 
S3=SHEAR(J-1) *L(J) *L(J) / (2*E*XI(J)) 
SLOPE (J) =S1-S2-S3 
Y1=SLOPE(J) *L(J) 
Y2=T(J-1) *SLOPE(J-1) *(L(J) **3) /(3*E*XI(J)) 
Y3=MOM(J-1) *(L(J) **2) /(2*E*XI (J) ) 
Y4=SHEAR(J-1) *(L(J) **3) /(3*E*XI(J)) 
Y (J) =Y¥ (J-1) -¥1+Y¥2-Y3-Y4 
SHEAR (J) =SHEAR(J-1) +( (WT (J) *WSQ*Y(J)) /(32.174*12) ) 
MOM1=SHEAR(J-1) *L(J) 
MOM2=T (J-1) *(¥(J-1)-Y(J)) 
MOM (J) =MOM(J-1) +MOM1-MOM2 

030 CONTINUE 
IF(I.EQ.2) GOTO 200 
BS=SHEAR (17) 
BM=MOM (17) 
BSLOPE=SLOPE(17) 
BY=Y (17) 
I=2 

C 

CER KKK KKK KEKE KKK KKK KKK KKK KKK KKK KK KKK KKA KK KKK KKK 


C NEW INITIAL CONDITIONS FOR SECOND COEFFICIENTS 


SLOPE(1)=0 
Y(1)=1 
Cc 
Crk a a KKK aKa Ki te a KKK KK KKK KKK KKK KKHMKE KK KKK 
GOTO 020 
200 AS=SHEAR(17) 
AM=MOM(17) 
ASLOPE=SLOPE(17) 
AY=Y(17) 
Cc 
C DELTAl VARIABLE IS FOR THE ARTICULATED ROTOR 
C DELTA2 IS FOR THE HINGELESS ROTOR 
Cc 
C WHERE DELTA GOES TO ZERO IS THE NATURAL FREQUENCY 
Cc 
DELTA 1=(AM*BY) - (BM*AY) 
DELTA2=(ASLOPE*BY) -(BSLOPE*AY) 
c 
50 FORMAT (1X,1PE13.6,1X,1PE12.5,1X,1PE12.5,1X, 1PE12.5,1X,1PE12.5, 


$1X,1PE12.5) 


78 


51 


400 


FORMAT (1X,F10.3,2X,F10.3) 
WRITE(*,50) W,DELTA1,AM, AY, BM, BY 
WRITE(35,50) W,DELTA1,AM, AY, BM, BY 
CONTINUE 


STOP 
END 


iS 


APPENDIX B: 


First Natural 


Operating rpm 


Frequency 
(rad/s) 


-071522E+02 
.071522E+02 
-071522E+02 
-071522E+02 
.071522E+02 
~071522E+02 
-071522E+02 
-071522E+02 
-071522Et02 
-O71522E+t+02 
-071522E+02 
-O71522Et+02 


ee) 


Frequency at 
107.1522 


(1022 rpm) 


Frequency 


759 (79.4118 


Delta 


6.82392E-01 
5.74431E-01 
4.66470E-01 
3.58510E-01 
2.50549E-01 
1.42588E-01 
3.46273E-02 
=7.33335E-02 
-1.81294E-01 
-2.89255E-01 
-3.97216E-01 
=9.05177 E-01 


which Delta 


0.0 


MYKLESTAD NATURAL FREQUENCIES 


ra 


Am 


bat ete ttt et 


d/s) 


-79121E+04 
-79121E+04 
-79121E+04 
-79121E+04 
-79121E+04 
-79121E+04 
~-79121E+04 
-79122E+04 
~-79122E+04 
-79122E+04 
-79122E+04 
-79122E+04 


-79121E+04 


80 


Coefficients 
Ay Bm 
7.48559E+00 -1.12384E+06 
7.48560Et+00 -1.12384E+t06 
7.48560E+00 -1.12384E+06 
7.48560E+00 -1.12384E+06 
7.48560E+00 -1.12384E+06 
7.48560E+00 -1.12384E+06 
7.48560EtO0O -1.12384E+06 
7.48560E+00 -1.12384E+06 
7.48560E+00 -1.12384E+06 
7.48560E+00 -1.12384E+06 
7.48560E+00 -1.12384E+06 
7.48561E+00 -1.12384E+06 
7.48560 -1.12384E+06 


By 


-4.69659E+02 
-4.69659E+02 
-4.69659E+02 
—-4.69659E+02 
-4.69659E+02 
-4.69659Et+02 
-4.69659E+02 
-4.69659E+02 
-4.69659E+02 
-4.69659E+02 
-4.69659E+02 
—4.69659E+02 


-4.69659E+02 


2nd Natural Frequency 


Operating rpm 759 (79.4118 


Frequency 
(rad/s) 


3.381635E+02 
3.381635E+02 
3.381636Et02 
3.381637E+02 
3.381638E+02 
3.381639E+02 
3.381639E+02 
3.381640E+02 
3.381641E+02 
3.381642Et02 
3.381643E+02 
3.381643E+02 
3.381644E+02 
3.381645E+02 
3.381646E+02 
3.381647E+02 
3.381647E+02 
3.381648E+02 
3.381649E+02 


Delta 


i 
WAOHKkY He ew 


t 
N= =e OD Wm 


-15818E+01 
-90265E+01 
-64711E+01 
-39158E+01 
-13605E+01 
-80514E+00 
-24980E+00 
-69446E+00 
-13912E+00 
-41623E+00 
-97157E+00 
-52692E+00 
.08227E+00 
-16376E+01 
-41930E+01 
-67483E+01 
~-93037Et+01 
-18591E+01 
2. 


44144E+01 


Frequency at whci Delta = 


338.1641 


(3229 rpm) 


OG. 


0 


rad/s) 


Am 


3.50008E+05 
3.50008E+05 
-SOO008E+05 
- 50009E+05 
- 50009E+05 
- 50009E+05 
- 50010E+05 
-50010E+05 
-SO01OEt+O5 
-50010E+05 
-90011E+05 
-50011E+05 
-50011E+05 
-50011E+05 
-50012E+05 
-50012E+05 
-50012E+08 
-50012E+05 
-50013E+05 


WQWWWWW WWW www Ww 


0 


3.50010E+05 


Coefficients 

Ay Bm 
8.05647E+01 -5.16530E+06 
8.05647E+01 -5.16530E+06 
8.05647E+01 -5.16530E+06 
8.O05648E+01 -5.16531E+06 
8.05648E+01 -5.16531Et06 
8.05649E+01 -5.16531E+06 
8.05649E+01]1 -5.16532E+06 
8.05650E+01]1 -5.16532E+06 
8.O05650E+01 -5.16532Et+06 
8.05651E+01 -5.16532E+06 
8.05651E+01 -5.16533E+06 
8.05652E+0)1 -5.16533E+06 
8.05652E+01 -5.16533E+06 
8.05652E+01 -5.16534E+06 
8.05653E+01 -5.16534E+06 
8.05653E+01 -5.16534E+06 
8.05654E+01 -5.16534E+06 
8.05654E+01 -5.16535E+06 
8.05655E+01 -5.16535E+06 
8.05650E+01]1 -5.16532E+06 


81 


By 


=; 
=1", 
she 
=i. 
-l. 
gh er 
oe 
= 
= 
=I. 
al Ni 
aoa 
-1. 
=i, 
-l. 
=i; 
=i. 
=4". 
=1. 


18895E+03 
18895E+03 
18895E+03 
18895Et03 
18895E+03 
18895E+03 
18895E+03 
18895E+03 
18895E+03 
18895E+03 
18895E+03 
18895E+03 
18895E+03 
18895E+03 
18895E+03 
18895E+03 
18895E+03 
18895E+03 
18895E+03 


- 18895E+03 


3rd Natural Frequency 


Operating rpm 759 (79.4118 


Frequency 
(rad/s) 


a i i i Mie i Be i i i | 


-583725EtO02 
-583725Et02 
~-583725E+02 
~-583725E+02 
~-983725Et02 
-o88 7256102 
-583725E+O02 
-583725E+02 
~-583725Et02 
~-583725E+02 
~-583725Et+O2 
-583725E+02 
-583726EtO02 
-583726E+02 
/ 983 726Et02 
~583726Et02 
-583726Et+02 
-583726Et+02 


Frequency at 


758.3725 


(7242 rpm) 


Delta 


m= WN) Oe BP Re vv 


.20479E+01 
.O0033EtO1 
-79588E+0O1 
-59143E+01 
. 38698Et+O1 
.18253E+01 
.78079E+00 
.73628E+00 
-69177Et00 
.64725E+00 
-60272E+00 
-~4. 
-2. 
“4. 
=6.; 
-8. 
alee 
=i. 


41768E-01 
48627E+00 
53081E+00 
57532E+00 
61982E+00 
O6643E+01 
27089E+01 


which Delta 


0. 


0 


rad/s) 


Am 


8.24037E+06 
8.24038E+06 
8.24038 E+ 06 
8. 24038E+06 
8.24038E+06 
8.24038E+06 
8.24038E+06 
8.24038E+06 
8.24038E+06 
8.24038Et+06 
8.24038Et+06 
8.24038E+06 
8.24038 E+06 
8.24038E+06 
8.24038E+06 
8.24038Et+06 
8.24038E+06 
8.24038E+06 


0 


8.24038E+06 


Coefficients 

Ay Bm 

8.12285E+O2 -7.04041E+t+07 
8.12285E+tO2 -7.04041E+07 
8.12286E+02 -7.04041Et07 
8.12286E+O02 -7.04041E+07 
8.12286E+O02 -7.04041E+07 
8.12286E+O2 -7.04041E+t+O?7 
8.12286EtO2 -7.04041E+t+07 
8.12286E+O2 -7.04041Et0O7 
8.12286E+tO2 -7.04041Et07 
8.12286EtO2 -7.04041E+07 
8.12286E+O2 -7.04041Et+07 
8.12286EtO2 -7.04041Et07 
8.12286E+O2 -7.04041E+07 
8.12286EtO2 -7.04041E+07 
8.12286EtO2 -7.04041E+07 
8.12286EtO2 -7.04041E+07 
8.12286Et+O2 -7.04041E+07 
8.12286E+O02 -7.04041Et+0O7 
8.12286EtO2 -7.04041E+07 


er 


By 


=6. 
-6. 
=65 
=6. 
=67 
= Gre 
-6. 
-6. 
me) o 
=6 
-6. 
-6. 
=62 
=6. 
-6. 
=6. 
=6 
=6. 


-6. 


94000Et+03 
94000Et+03 
94000Et+03 
94000E+03 
94000E+03 
94000E+03 
94000E+03 
94000Et03 
94000E+03 


-94000E+03 


94000Et+03 
94000Et+03 
94000Et03 
94000E+03 
94000E+03 
94000Et03 


-94000E+03 


94000E+03 


94000E+03 


4th Natural Frequency 


Operating rpm 759 (79.4118 


Frequency 
(rad/s) 


a a ee ee 


-448090E+03 
-448091E+03 
-448092E+03 
~-448092E+03 
-448093E+03 
-448094E+03 
~-448095E+03 
-448096E+03 
-448096E+03 
-448097E+03 
-448098E+03 
~-448099E+03 
~-448100E+03 
-448100E+03 
1. 
1. 
1. 
1. 


448101E+03 
448102E+03 
448103E+03 
448104E+03 


Frequency at 


1448.100 


(13828 rpm) 


-4 
=o 
=3. 
=32., 
-2. 
me 
-2. 
=1; 
-1. 
=1.. 
=6:. 
=o 5 
=2. 


Delta 


~-21455E+04 
88022E+04 
54589E+04 
21156E+04 
B87722E+04 
-54288E+04 
20854E+04 
87420E+04 
53985E+04 
20551E+04 
71160E+03 
36809E+03 
02456E+03 
1.31898E+03 
4.66255E+03 
8.00613E+03 
1.13497E+04 
1.46934E+04 


which Delta 


0.0 


rad/s) 


NONNNNNNNN NN NN NN AY YH AY 


Am 


-77806E+08 
- 77807E+08 
- 77808E+08 
-77809E+08 
-77810E+08 


e 


77811E+08 


- 77812E+08 
-77813E+08 
-77814E+08 
-77815E+08 
-77816E+08 
-77817E+08 
-77818E+08 
-77819E+08 
- 77820E+08 
~-77821E+08 
~77822E+08 
»-77822E+08 


-77819E+08 


83 


Coefficients 


Ay 


1.18496E+04 
1.18497E+04 
1.18497E+04 
1.18497E+04 
1.18498E+04 
1.18498E+04 
1.18498E+04 
1.18499E+04 
1.18499E+04 
1.18499E+04 
1.18499E+04 
1.18500E+04 
1.18500E+04 
1.18500E+04 
1.18501E+04 
1.18501E+04 
1.18501E+04 
1.18502E+04 


1.18500E+04 


Br 


=1. 752506709 
=]. 7525158109 
-1.75251E+t09 
=1. 75252E+09 
-1.75253E+09 
-1.75253E+09 
-1.75254E+09 
-~1.75254E+09 
=). 752556709 
=i /5259E8409 
-~1.75256E+09 
=1.75257E109 
-1.75257E+09 
-1.75258E+09 
-1.75258E+09 
=) 752599E+09 
-1.75259E+09 
-1.75260E+09 


-1.75258E+09 


By 


-7 
-7 
-7 
-7 
-7 
-7 
-7 
-7 
-7 
-7 
-7 
-7 
-7 
-7 
-7 
-7 
-7 
-7 


~-47518E+04 
»-47520E+04 
-47521E+04 
-47523E+04 
-47525E+04 
-47527E+04 
»-47529E+04 
»-47531E+04 
~-47532E+04 
~-47534Et04 
»-47536E+04 
-47538E+04 
-47540E+04 
~-47542E+04 
~-47544E+04 
-47545E+04 
»-47547E+04 
»-47549E+04 


-47541E+04 


First Natural 


Operating rpm 


Frequency 
(rad/s) 


1 WO WO WO 1 0 WO WO 1 WO WO WO 10 


-987550E+01 
-987551E+t0O1 
~987552Et+01 
-987553Etr01 
-987553E+01 
-987554E+01 
-987555Et01 
-987556E+01 
-987557E+01 
-987557E+01 
-987558E+01 
790 7999E 401 
-987560E+01 


Frequency at 


99°57 87 S57 


(954 rpm) 


Frequency 


(Light Blade) 


759 (79.4118 rad/s) 


Delta 


NJ BP ND RD | WwW 


-75536E-01 
- 32246E-01 
.88956E-01 
-45666E-01 
.02376E-01 
-59087E-01 
-15797E-01 
-25066E-02 
-92166E-02 
eh. 
-73633E-02 
=e 
-l]. 


40733E-02 


00653E-01 
43943E-01 


which Delta 


0. 


0 


Am 


DANAANAAAAAAA NH 


-90470Et03 
.90470E+03 
.90470Et03 
-90470E+03 
-90470E+03 
-90471E+03 
-90471E+03 
-90471Et03 
-90471E+03 
-90471E+03 
-90471E+03 
-90471E+03 
-90471E+03 


-90471Et03 


Coefficients 

Ay Bm 
2.51694E+00 -4. 
2.51694E+00 -4. 
2.51694E+00 -4. 
2.51694Et+00 -4. 
2.51694E+00 -4. 
2.51694E+00 -4. 
2.51694E+00 -4. 
2.51694E+00 -4. 
2.51694E+00 -4. 
2.51695E+00 -4. 
2.51695E+00 ~4. 
2.51695E+00 -4. 
2.51695E+00 ~-4. 
2.51695 =-4., 


84 


13426Et05 
13426E+05 
13426E+05 
13426E+05 
13426E+05 
13426E+05 
13426E+05 
13426E+05 
13426E+05 
13426E+05 
13426E+05 
13426E+05 
13426E+05 


13426E+05 


By 


-1.50704E+02 
-1.50704Et02 
-1.50704E+02 
-1.50704E+02 
-1.50704E+02 
-1.50704E+02 
-~1.50704E+02 
-1.50704E+02 
-1.50704E+02 
-1.50704E+02 
-1.50704E+02 
-1.50704E+02 
-1.50704E+02 


-~1.50504E+02 


2nd Natural Frequency (Light Blade) 


Operating rpm 759 (79.4118 rad/s) 


Frequency 
(rad/s) 


abr b hb bhbhbhbbhbbhdhbbhbhbhbhbh bb b&b & 


-004961E+02 
-004961Et0O2 
.004961E+02 
-004961Et02 
.004961Et02 
-CO4962Et02 
-004962Et02 
-004962E+02 
-004962Et+02 
-004962Et02 
-004962E+02 
-004962Et02 
-004962Et+02 
-004962E+02 
-004962E+02 
-004962Et02 
-004962E+02 
-004962Et+02 


Delta 


(elie et tt 
WO > = WwW & OV «)} © = 


-03982E+t00 
~-98344E-01 
-56864E-01 
-15385E-01 
-73906E-01 
324 268-01 
-90947E-01 
~-94674E-02 
~20121E-02 


2.33492E-01 


Rm eH OH WwW 


~74971E-01 
-16451E-01 
~57930E-01 
-99410E-01 
-40889E-01 
-08237E+00 
-22385E+00 
-36533E+00 


Frequency at which Delta 


400.4962 


(3824 rpm) 


QO. 


0 


Am 


ee ee ee) 


-91468EtO05 
~91468E+05 
-91468E+05 
-91468E+05 
-91468E+05 
-91468E+05 
-91468E+05 
~-91468E+t05 
-91468E+05 
-91468E+05 
-91468E+05 
~-91468Et05 
-91468E+05 
-91468E+05 
-91468E+05 
-91468E+05 
-91468Et05 
~-91468E+05 


-91468E+05 


Coefficients 


Ay 


3. OSs 55ET0)] 
3.05355Et01 
3.05355E+01 
3.05355Et01 
3.05355Et01 
3.05355Et01 
3.05355E401 
3.05355E+01 
3.05356Et0] 
3.05356E+t01 
3.05356E+01 
3.05356E+01 
3.05356Et0O1 
3.05356Et01 
3.05356E+t01 
3.05356E+01 
3.05356E+01 
3.05356Et0O1 


3.05356Et0O] 


85 


Bm 


- 


-2.91981E+06 
-2.91981E+06 
-2.91981Et06 
-2.91981E+06 
-2.91981E+06 
-2.91981E+t06 
-2.91981Et06 
~2.91981E+06 
-2.91981E+06 
-2.91981E+t06 
-2.91981E+06 
-2.91981E+t06 
-2.91981Et06 
-2.91981E+06 
-2.91982E+06 
-2.91982E+06 
-2.91982Et06 
=2.,91982E106 


-2.91981E+06 


By 


-4 
-4 
-4 
-4 
-4 
-4 
-4 
-4 
-4 
-4 
-4 
-4 
-4 
-4 
-4 
-4 
-4 
-4 


-65656E+02 
-65656Et02 
-65656E+02 
-65656Et02 
-65656E+02 
-65656E+02 
-65656E+02 
-65656E+02 
-65656Et02 
-65656E+02 
-65656Et02 
-65656E+02 
-65656E+02 
-65656Et+02 
-65656E+02 
-65656E+02 
-65656E+02 
.65656E+02 
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3rd Natural Frequency 


(Light Blade) 


Operating rpm 759 (79.4118 rad/s) 


Frequency 
(rad/s) 


1.062360E+03 
1.062361E+03 
1.062362Et03 
1.062362E+03 
1.062363E+03 
1.062364E+03 
1.062365E+03 
1.062365Et03 
1.062366E+03 
1.062367E+03 
1.062368E+03 
1.062369E+03 
1.062369E+03 
1.062370E+03 
1.062371E+03 


Frequency at 
1062.366 


(10145 rpm) 


Delta 


-49707E+03 
-31505E+03 
-13303E+03 
-51010E+02 
-68988E+02 
-86965Et02 
-04941E+02 
-22916E+02 
-08897E+01 
-~1.41137E+02 
-3.23166Et02 
-5.05195E+02 
=6.6/2255102 
-8.69256E+02 
-1.05129E+03 


&M 2 UW OH ee 


which Delta 


0.0 


SSN N™NWNNNNNN SYN SS 


-06185E+06 
-06187E+06 
-06189E+06 
-06192E+06 
-06194E+06 
-06197E+06 
-06199E+06 
-06201E+06 
-06204E+06 
-06206E+06 
-06209E+06 
-06211E+06 
-06213E+06 
-06216E+06 
-06218E+06 


-06204E+06 


86 


Coefficients 


Ay 


4.97993E+02 
4.97995E+02 
4.97996E+02 
4.97997E+02 
4.97999E+02 
4.98000Et+02 
4.98001E+02 
4.98003E+02 
4.98004E+02 
4.98005E+02 
4.98007E+02 
4.98008Et02 
4.98009E+02 
4.98011E+02 
4.98012E+t+02 


4.98004Et02 


” 


Bm 


-6.33427E+07 
-6.33429E+07 
-6.33431E+07 
-6.33433E+07 
-6.33435E+07 
-6.33437E+07 
-6.33439E+07 
-6.33441E+07 
-6.33443E+07 
-6.33444E+07 
-6.33446E+07 
-6.33448E+07 
-6.33450E+07 
-6.33452E+07 
-6.33454E+07 


-6.33443E+07 


By 


-4.46686E+t+03 
-4.46687E+03 
-4.46688E+03 
-4.46689E+03 
-4.46690E+03 
-4.46691E+03 
~4.46692E+03 
-4.46693E+03 
-4.46694E+03 
-4.46695E+t03 
-4.46696E+t03 
-4.46697E+03 
-4.46698E+03 
-4.46699E+t03 
-4.46700E+03 


-4.46694E+03 


4th Natural Frequency 


Operating rpm 759 


Frequency 
(rad/s) 


-O90913Et0O3 
-O90913E+03 
-090914E+03 
.090915E+03 
.090916E+03 
.090917Et03 
.090917E+03 
.090918E+03 
.090919E+03 


NNNNN NN ND PY 


Frequency at 
2090.917 


(19967 rpm) 


Delta 


me 
-15527E+04 


=1 


-8. 
Sor 
=1. 
-29533Et+03 


1 


4. 
7. 
iz 


47646E+04 


34069E+03 
12869E+03 
91669E+03 


50736E+03 
71941E+03 
09315E+04 


which Delta 


0. 


0 


Am 


NNNN NN NY ND LH 


=) 


(Light Blade) 


(79.4118 rad/s) 


~-52950E+08 
-52951E+08 
-52952E+08 
-52952E+08 
-52953E+08 
-52953E+08 
-52954E+t+08 
-52955E+08 
-52955E+08 


~-52953E+08 


a7 


Coefficients 


Ay 


9.15755E+03 
9.15757E+03 
9.15759E+03 
9.15761E+03 
9.15762E+03 
9.15764E+03 
9.15766E+03 
9.15768E+03 
9.15769E+03 


9.15764E+03 


Bm 


-1.68095E+09 
-1.68096E+09 
-1.68096E+09 
-1.68097E+09 
-1.68097Et09 
-1.68097Et09 
-1.68098E+09 
-1.68098E+09 
~1.68098E+09 


-1.68097E+09 


By 


=6% 
=6. 
al 
=G: 
=64 
= 6:4 
~6. 
~6. 
=67. 


08555E+04 
08556E+04 
08557E+04 
O08558Et04 
08559E+04 
O08560E+04 
0O8561E+04 
08563E+04 
08564E+04 


-60.8560E+04 


APPENDIX CC: MODE SHAPE FORTRAN CODE 


xakKaKSK 


x 
x 
& 


This Program is used to calculate the mode shapes once 
the natural frequencies are determined by the Myklestad 
determinant method 


ake 


CUO Gl G1 Ola Oe il ee @. 


IMPLICIT DOUBLE PRECISION (A-H,0-2Z) 

DOUBLE PRECISION RSTA(24) ,DMASS(24),EI(24) ,CENT(24) , SHEAR(24) 
DOUBLE PRECISION DMOM(24),SLP(24),DEFL(24) ,ix(24),am(24), 
4dx(24),mod,pi,p 

OPEN (101, FILE='MODE. OUT’ , STATUS=/ UNKNOWN’ ) 


~ 


pi = 3.141592654 ! PI 

gq = 32.174 ! GRAVITATIONAL CONSTANT 

e = 0.133 ! HINGE OFFSET RATIO 

RAD = 68.0 ! ROTOR RADIUS 

RO = e*RAD ! HINGE POINT 

C = 6.5 ! CHORD 

to = .13 ! ROOT SECTION THICKNESS (t/c) 

tt = .13 ! TIP SECTION THICKNESS (t/c) 
= 0) ! DENSITY OF SPRUCE lbm/in‘*3 

mod = 1305000.0 ! MODULUS OF ELASTICITY SPRUCE 


PRINT *, ‘ENTER ROTOR RPM AND NATURAL FREQUENCY ’ 
READ *,rpm,wrpm 

RV=rpm*2.0%*pi/60.0 

W=wrpm*2.0*pi/60 


THE FOLLOWING ARRAYS ARE NUMBERED FROM BLADE TIP TO ROOT; 


SEGMENT RADIUS (IN) 
DATA RSTA/65.6,62.0,59.6,57.2,54.8,52.4,50.0,47.6,45.2,42.8, 
$40.4,38.0,35.6,33.2,30.8,28.4,26.0,23.6,21-2,18.8, 1684, lane 
$11.6,9.2/ 


SEGMENT SPACING (IN) 
DATA adx/4.8,23*2.4/ 


SEGMENT WEIGHT (LBS) 
DATA dm/.353,20%.177,.302,.252,.470/ 


SEGMENT FLAPWISE MOMENT OF INERTIA, IXX (IN*4) 
DATA ix/21*.478,.48, .400, .845/ 


DO 10 J=24,1,-1 
x = (J-1) *dx(LL) +0.5*dx (LL) 
RSTA(25-J)= x 
t = ((RAD-RO-x) *(t0O-tt) /(RAD-RO) +tt) *C 
QGm = (2.0%6.60*0.2850+2.0*t*0.1250) *p*dx (LL) 
DMASS(25-J) = dm/g 
i = (2.0*6.85%.2850) *((.5*(t-.1425) ) **2.0) 
EI(25-J) = i*mod 


10 CONTINUE 


CENTRIFUGAL FORCES COMPUTED ALONG THE BLADE 
DMASS(1)=dm(1)/g 


CENT (1) =DMASS(1) * (RO+RSTA(1)) *(RV#**2) /12 
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30 


! 


30 


40 


60 


DO 30 J=2,24 

DMASS(J)=dm(J)/G 

EI (J)=ix(J) *mod 

CENT (J) =CENT (J-1) +DMASS (J) * (RO+RSTA(J) ) *(RV**2) /12 
CONTINUE 


ITERATION LOOP 


SLP(1) = 0.0 ! BC’S AT TIP FOR Ay AND Am 
DEFL(1) = 1.0 
DMOM(1) = 0.0 
SHEAR(1) = 0.0 
DO 40 KK=1,2 
DO 50 LL=2,24 ! MYKLESTAD INTEGRATION 


SLP1=SLP(LL-1) * (1+CENT(LL-1) * (dx (LL) *#*2) /2/EI (LL) ) 
SLP2=DMOM(LL-1) *dx (LL) /EI (LL) 
SLP3=SHEAR (LL-1) * (dx (LL) **2) /2/EI (LL) 
SLP (LL) =SLP1-SLP2-SLP3 
DEFL1=SLP (LL) *dx (LL) 
DEFL2=CENT(LL-1) #SLP(LL-1) * (dx (LL) **3) /3/EI (LL) 
DEFL3=DMOM(LL-1) * (dx (LL) **2) /2/EI (LL) 
DEFL4=SHEAR(LL-1) * (dx (LL) **3) /3/EI (LL) 
DEF (LL) =DEFL(LL-1) -DEFL1+DEFL2-DEFL3-DEFL4 
SHEAR (LL) =SHEAR (LL-1) +DMASS (LL) * (w**2) *DEFL(LL) /12 
DMOM1=SHEAR(LL-1) *dx (LL) 
DMOM2=CENT (LL-1) * (DEFL(LL-1) -DEFL(LL) ) 
DMOM (LL) =DMOM (LL-1) +DMOM1-DMOM2 

CONTINUE 

IF (KK.EQ.1) THEN ! COEF. COMPUTATION 
Ay=DEFL(24) 
Am=DMOM (24) 

ELSE 
By=DEFL(24) 
Bm=DMOM (24) 


END IF 
SLP(1)=1.0 ! BC’S FOR By AND Ba 
DEFL(1)=0.0 

CONTINUE ! LOOP FOR By AND Bm 


DET=Am*By-Ay*Bm 
PRINT *,DET 


SLP (1) = ~Ay/By 
DEFL(1) = 1.0 
DMOM(1) = 0.0 
SHEAR(1) = 0.0 
DO 60 LL=2,24 ! MYKLESTAD INTEGRATION 


SLP1=SLP (LL-1) *(1+CENT(LL-1) * (dx (LL) **#2) /2/EI (LL) ) 
SLP2=DMOM(LL-1) *dx (LL) /EI (LL) 
SLP3=SHEAR(LL-1) * (dx (LL) #*2) /2/EI (LL) 
SLP (LL) =SLP1-SLP2-SLP3 
DEFL1=SLP (LL) #dx (LL) 
DEFL2=CENT(LL-1) *SLP(LL-1) * (dx (LL) **3) /3/EI (LL) 
DEFL3=DMOM (LL-1) * (dx (LL) *#2) /2/EI (LL) 
DEFL4=SHEAR (LL-1) * (dx (LL) #*3) /3/EI (LL) 
DEFL (LL) =DEFL(LL-1) -DEFL1+DEFL2-DEFL3-DEFL4 
SHEAR (LL) =SHEAR (LL-1) +DMASS (LL) * (w**2) *DEFL(LL) /12 
DMOM1=SHEAR (LL-1) *dx (LL) 
DMOM2=CENT (LL-1) * (DEFL(LL-1) -DEFL(LL) ) 
DMOM (LL) =DMOM ( LL- 1) +DMOM1-DMOM2 

CONTINUE 
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DO 100 J=1,24 
WRITE (101,*)J,RSTA(25-J) ,DEFL(25-J) 

100 CONTINUE 
kkkkkkkhk kkk khkkkkkkhkhhhkhkhkhkhkhkhkhkkkakkhkkkkhkkhkhkahkhhkkhkhkhhkhkakhhkkkkkhake 
k 
* This program was written by Lt Matt Avila 
* Modifications were made by Lt J.L. Vandiver 
k 
Eee RPeESCESESE ER SSEESE SE EERE RES ESE SSE REESE REESE EES EPEC EPEC Eee ee ee ee 


END 
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tHelk opter Design 


PRELIMINARY POWER CALCULATIONS 
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APPENDIX E: BLADE TIP PITCH ANGLE WORKSHEET 
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Pitch Angle Determination. 
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